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' We propose a trick which enables one to incorporate the Gribov horizon 

. into the Schwinger-Dyson equation in Landau and Coulomb gauge Yang-Mills 

Qh. theory, using the Gribov- Zwanziger framework with the horizon term. We find 

I a family of solutions parameterized by one-parameter wji{0) which was over- 

looked so far by assuming to be zero implicitly. The family includes both the 
. scaling and decoupling solutions, and specification of the parameter discrim- 

inates between them. In the Landau gauge we discuss a possible decoupling 
solution satisfying the Kugo-Ojima criterion for color confinement. 
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1 Introduction 



The Green functions are quantities of fundamental importance in quantum field the- 
ories. In quantum Yang-Mills theory [1], the Green functions are gauge-dependent 
quantities and defined only when the gauge fixing condition is imposed, as far as the 
local Green functions are concerned. In other words, the result of calculating the 
Green functions can be affected by details of the explicit procedure for gauge fixing, 
in sharp contrast to the gauge-invariant quantities which do not depend on the gauge 
fixing adopted and should take the same value even if they are calculated in different 
gauges. However, the gauge-invariant operators such as the Wilson loop operator and 
the Polyakov line operator are necessarily non-local. 

The quantum Yang-Mills theory in the Landau gauge ds^ = as a manifestly 
CO variant gauge is naively defined by the functional integral: 



y" \d£^\b{ds^) det(-9DK])) exp{-SyMK]}, (1-1) 



where Sym is the Yang- Mills action and —dD\s^\ is the Faddeev- Popov operator. 
However, the Landau gauge fixing ds^ = can not fix the gauge uniquely. This is 
because each gauge orbit intersects the gauge fixing hypersurface F := ; ds/ = 0} 
many times. The unique representative from each gauge orbit can not be chosen by 
imposing = alone. There are many representatives called the Gribov copies [2]. 

In order to avoid the Gribov copies, Gribov proposed to restrict the functional 
integral to the 1st Gribov region Q: 

^Gribov := Jjd^]S{ds^) det{-dD[j^])) exp{-^yM[^]}, (1.2) 

where the Gribov region is defined by 

n ■= {j^; ds^ = 0, -dD[£^] > 0} C F. (1.3) 

The 1st Gribov region is a bounded and convex region including the origin {s^ = 0}. 
In fact, -d^Df,[^ = 0] = -d^,d^ > 0, i.e., {s/ = 0} e Q. The boundary of Q is 
called the Gribov horizon: 

dn := dj^ = 0, -dD[£/] = 0}. (1.4) 

He predicted that the resulting Green functions exhibit unexpected behavior in the 
deep infrared (IR) region and that they play the essential role in confinement. 
We define the gluon 2-point function (full or complete propagator) by 



and the ghost propagator by 



^2 / 



(« = 0), (1.5) 



G^\k) := -5^^^, (1.6) 

where the free case corresponds to F{k^) = 1 and G{k'^) = 1. Gribov predicted their 
behaviors in the deep IR region A;^ <^ 1: 

-fc^~(PFW^°' ^(^')- IT Too (eiO), (1.7) 
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where M is a constant with mass dimension called the Gribov mass. This power like 
behavior should be compared with the ultraviolet (UV) behavior with the logarithmic 
corrections. 

The Gribov prediction was investigated more elaborately by solving the coupled 
Schwinger-Dyson (SD) equation for the gluon and ghost propagators. The dressing 
functions F and G are characterized by the power behavior with exponents a and /3 
respectively: 

F{k^)=Ax{k^Y, G{k^) = Bx{ky, a + 2/3 = 0, 0<A,B<oo. (1.8) 

By the scaling relation a + 2/3 = 0, a single exponent k is enough for characterizing 
the IR behavior: 

a = 2K>l, (3 = -K<0, 1/2<K<1 (Gribov K = l). (1.9) 

The Gribov prediction corresponds to the limiting value k, = 1. This result, gluon 
suppression and ghost dominance, leads to the running coupling constant with a 
non-trivial IR fixed point: 

g^{k):=g^F{e)G^{e)^{0<)g^AB^{<oo) (P^O). (1.10) 

This solution is called the scaling solution. The gluon propagator F{k'^)/k'^ vanishes 
in the IR limit k^ | 0, while the ghost propagator becomes more singular than the free 
case in the IR region, or the ghost dressing function G{k^) diverges, i.e., 

G(0) = oo. (1.11) 

See the excellent review by Alkofer and von Smekal [3] for details. 

This IR behavior was considered to be reasonable from the viewpoint of color 
confinement. Due to Kugo and Ojima [1], all color non-singlet objects can not be 
observed or confined, in other words, only color singlet objects are observed, if a 
criterion u{0) = —1 is satisfied in the Lorentz covariant gauge (a sufficient condition 
for color confinement). It is shown in [5] that in the Landau gauge, the Kugo-Ojima 
criterion for color confinement u{0) = —1 is equivalent to the divergent ghost dressing 
function G{0) = oo , since in the Landau gauge 

G{0) = [1 + u{0)]~\ (1.12) 

In this paper, we point out that this relation is not exact and must be used with more 
care, since m(0) = —1 does not necessarily mean G{0) = oo. 

Until 2006, it seemed that the scaling solution has been confirmed by the SD 
equation, the functional renormalization group equation and numerical simulations 
on lattice. This lead to the ghost dominance picture for color confinement. 

So far so good. However, these results are questioned by the Orsay group at Uni- 
versite de Paris Sud. By careful analyses of the SD equation, socalled the decoupling 
solution was discovered [H]: 

F{e) =A' X {kY, G{k^) = B'x {kY, = 1, /?' = 0, < A', B' < oo, 

(1.13) 

^ Note that the Kugo-Ojima theory for color confinement does not take into account the Gribov 
problem and is based on the usual BRST formulation where the exact color symmetry and the 
well-defined BRST charge are assumed. 



2 



which lead to the running couphng going to zero in the IR hmit: 

g\k) := g^F{e)G\k^) ~ g^A'B'^k^ _^ 0), (1.14) 

although its possibility was mentioned also in [7j. Moreover, reexaminations of nu- 
merical simulations on large lattices [Hl[9l[10l[ll],[T2],[l3l[llj , functional renormalization 
group equation [15] and other methods [T8],[l9] seem to support this result. The de- 
coupling solution implies that the gluon propagator goes to the non-zero and finite 
constant in the IR limit, while the ghost propagator behaves like free, namely, the 
ghost dressing function G(0) is non-zero and finite in the IR limit: 

< G(0) < oo. (1.15) 

In the decoupling solution, the gluon decouples below its mass scale and the ghost 
is still dominant, although the ghost dominance in the decoupling solution is weaker 
than that in the scaling solution. 

In this paper, we take into account the existence of the Gribov horizon by making 
use of the Gribov- Zwanziger theory [20 11211 [22] : 

Zgz := j [d£^]5{d''£^^) det{K[£^]) exp{-5yA/[^] -i j d^xh{x)}, (1.16) 

where K is the Faddeev-Popov operator K[^] := —d^D^ = —d^{d^ + g^^x) and 
h{x) = h[s?/]{x) is the Zwanziger horizon function. Here the parameter 7 called 
the Gribov parameter is determined by solving a gap equation, commonly called the 
horizon condition: for D-dimensional Euclidean SU(N) Yang-Mills theory, 

{h{x))Gz = {N' - 1)D. (1.17) 

The horizon function plays the role of restricting the integration region inside the 
Gribov horizon. However, the exact form of the horizon function in this sense is 
not known and the choice of the horizon function is not unique at present. Some 
arguments on this point are given later. The first choice of the horizon function 

is [20] 

h{x) = I d^ygf^^'^^^^ix)iK-'f^ix,y)gf^^^^fiy). (1.18) 
The second choice is [22] 

h{x) = J d''yD,[s^f^{x){K-'f^{x,y)D,[j^f''{y). (1.19) 

In any case, inclusion of the horizon term makes the theory non-local. 

The SD equations for Green functions do not change their form even in the pres- 
ence of the Gribov horizon, since the integrand of the functional integration formula 
for Green functions vanishes at the Gribov horizon which is the boundary of the func- 
tional integration region due to the Faddeev-Popov operator. Hence, the solutions of 
the SD equation include both the solution with the Gribov horizon and the solution 
without restriction. Remarkably, it has been shown [T5|[T6] that the set of solutions of 
the coupled SD equation for gluon and ghost propagators is uniquely determined once 
a boundary value G(0) is given, corresponding to the scaling solution for G(0) = 00 
and the decoupling solution for < G(0) < 00. However, it is not yet examined how 
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these solutions are related to the Gribov horizon. Moreover, the horizon term does 
not uniquely fix the gauge, since there are still Gribov copies in the Gribov region. A 
one-parameter family of correlation functions are constructed in lattice gauge theory 
distinguished by a second gauge parameter B (Landau- -B gauge) [l7j. This uniquely 
specifies a representative from a gauge orbit and no further freedom in choosing a 
Gribov copy. In this paper, we give a trick to incorporate the horizon condition into 
the SD equation of the ghost propagator, which enables us to distinguish the solution 
associated with the particular choice of the horizon term. This is possibly used to dis- 
crimate the scaling and decoupling solutions. We show that both horizon terms allow 
the existence of one-parameter family of solutions parameterized by a real number 
wr{0) which has been assumed implicitly to be zero w_r(0) = in the previous in- 
vestigations. Therefore, wji{0) plays the role of an additional non-perturbative gauge 
parameter which uniquely specifies the solution. We consider both the unrenormal- 
ized and renormalized versions of the SD equation with the horizon condition being 
included. 

2 Schwinger-Dyson equation with the horizon con- 
dition being inserted 

In what follows, we consider the SU(N) Yang-Mills theory in dimensional Euclidean 
space. 

The Schwinger-Dyson (SD) equation for the ghost propagator {^^'i'^)^ in mo- 
mentum space is written in the following form. We follow the notation of [211I2S]- 



(2.1) 



which is obtained as the Fourier transform of 



= -{{d,D,[s^]^)\x)^\y)) + 5^^5^(x - y). 



(2.2) 



In the Gribov theory, this is derived from the identity: 




(2.3) 



while in the Gribov-Zwanziger theory, the same form is obtained from 




(2.4) 



where 



rtot 
'YM 



5'ym + 'S'gf+fp, 




(2.5) 
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Figure 1: Diagrammatic representation of {{g£^^ x 'io)^{g£^i, x ^)^)k, {{,9'^^ x 

X ^)^Yr\ {{gs^^ X <^)^(^?< X and x ^)^{gK x 

and the dot and the cross are defined as ■ S§ := s^^SS^ and x i^)^ := 
By using the relationship derived in 



^^{{g^, X = -^K^{k) = 6^^[u{e) + w{k% (2.6) 



the SD equation for the ghost dressing function G{k'^) is rewritten as 

G~\k'') = l+u{k^)+w{k''). (2.7) 

This identity was derived by [3] and also in [2S] based on a different method. Here 
two functions u and w; are defined from the modified 1-particle irreducible (ml PI) 
part as 



X^^^{k) := {{gs^, X ^)\gs^^ x 



5^^, (2. 



where u{k ) agrees with the Kugo-Ojima function usually defined by 



{{D,^)\g^, X ^f), := (^6,. - S^%{e). (2.9) 

The mlPI part is defined from the two-point function of the composite operators (See 
Fig. [ID 

{{g^, X ^^gK X ^f), = X^^^ik) + A^^^ik), (2.10) 

where 

K^{k) :=((^< X ^)V^)r(^^^"^),(^^(^< X (2.11) 

Following the idea of Gribov [2], we incorporate the horizon condition {h{0)) = 
[N"^ — 1)D into the SD equation. By substituting the horizon condition into the free 
(tree) part of the SD equation, the SD equation for the ghost dressing function reads 

G-\k^) = + <k^) + Hk% (2.12) 

which is equivalent to the SD equation for the ghost propagator: 

(^V^)- = -S^'k' J^^^l^ - ^r^y^'. (2.13) 
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2.1 The first horizon term 



The first horizon term fll.181) yields 

{h{0)) = - lim(((7< X ^)^{gs^, x 

= -{N^- 1) {Du{0) + w{0) - G{0)[u{0) + w{0)]^} , (2.14) 

and 

WO)) _ .„(o)_^+G-(0)-2 + G(0)^ p^^^, 



Then the SD equation reads 

^ G-(0)-^2 + G(0) _ ^ _ ^^^j ^ ^ p^^j 

Then, it is observed that the term — ^^(0) coming from the horizon condition exactly 
cancels the self-energy term u{k'^) at /c = in the right-hand side of the SD equation. 
In the deep IR limit k = 0, therefore, we have 

G-(0) = ^:M^l±^+f-l + lV-(0). (2.17) 



D \ D 

By solving this equation for G{0): 

G\0) - [2 - (D - l)w(0)]G(0) + 1-D = 0, (2.18) 

we have 

G(0) = 1-{D- l)u;(0)/2 + ^J[l -{D- 1)m;(0)/2]2 - 1 + D > 0, (2.19) 

and u{0) = -1 - w{0) + ^-^(0), i.e., 

1 



n(0) = -1 - 'u;(0) - -{2 - 3w{0) - ^12 + [2 - 3w{0)]'^}. (2.20) 

This implies that the horizon condition determines the boundary value G{0) in the 
ghost SD equation. Consequently, we have one-parameter family of solutions param- 
eterized by w{0). 

We consider G(0) and m(0) as functions of w{0). Both G(0) and u{0) are monoton- 
ically decreasing functions in w(0); G(0),m(0) ^ oo as w{0) —oo, while G{0) 
and n(0) —5/3 as w{0) +oo. For D = 4, in particular, G{0) = 3 and 
u{0) = -2/3 at w{0) = 0; G(0) = 2 and m(0) = -1 at w{0) = 1/2. 

Thus the scaling solution G(0) = cxo is obtained only when w{0) = — oo. Otherwise 
w(0) > — oo, the decoupling solution < G(0) < oo is obtained. Using a special 
value as an additional input w{0) = by an assumption [ni[21] or by an independent 
argument [221 G{0) is determined selfconsistently by solving the above SD equation 
as [2a[2S] 

G(0) = 1 + /D>0, m(0) = (-D±/D)/(L'- 1). (2.21) 

In particular, for D = A, 

G{0) = 3 > 0, u{0) = -2/3 {D = 4). (2.22) 

In order to obtain the scaling solution, the constant terms must cancel exactly or 
disappear at the k = limit on the right-hand side of the SD equation. This is what 
implicitly assumed, but not stated explicitly, as pointed out by [H]. 
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2.2 The second horizon term 

The second horizon function fll.l9p leads to 



(hiO)) = - limiiD^^nD,^)^), = -{N' - 1) {{D - l)n(O) - 1} , (2.23) 

and 

Then the SD equation ^"^(A;^) = 1 + n(A;^) + w{k'^) is rewritten as 

G-\k'') = ^[1 + m(0)] - m(0) + u{e) + u;(A:2). (2.25) 

In the deep IR hmit, the cancellation occurs for u(0) in the right-hand side of this 
equation: 

G'^'(0) = ;^[l + ^(0)] + MO), (2.26) 
which is rewritten in favor of G'(O) using G(0)~^ = 1 + m(0) + w(0) as 

G-i(O) = ^G-\Q) - ^w{0) + w{0). (2.27) 
This is solved for Z) 7^ 1 to give 

G{0) = w~\0), (2.28) 

and 

n(0) = -1. (2.29) 

Thus, we have one-parameter family of solutions parameterized by w{0). The scaling 
solution G{0) = cx) is obtained only when w{0) = 0. Otherwise < w{0) < 00, the 
decoupling solution 00 > G(0) > is obtained. It should be remarked that the Kugo- 
Ojima condition m(0) = —1 is always satisfied. El However, this does not immediately 
mean the enhancement of the ghost propagator, contrary to the usual claim found in 
literatures. 

If we require that two horizon conditions (12.141) and fl2.23p give the same result, 
then the relation 

u{0) = -l/2-w{0) (2.30) 

must be satisfied for any D. This implies that the common solution is found for any 
D 

G(0)=2, m(0) = -1, w{0) = -. (2.31) 



^ The fact that the horizon condition using this definition (|2.23p is equivalent to the Kugo-Ojima 
criterion (|2.29p has already been pointed out and it was checked to what extent the horizon condition 
holds in the numerical simulation in [30j . 
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3 Removing ultraviolet divergence and renormal- 
ization 

We have calculated the first horizon condition using 

{h{x)) = - lim(((7=< X X <^)^)fc. (3.1) 

However, the composite operator {gsz/^ x ^) is not multiplicative renormalizable due 
to operator mixing: 

gs^^x^ = Zc^'\gs^^ x + Zc^'\l - Zc)d^^R, (3.2) 

where ^ = Z^^'rfji. This implies that in the first choice one can not use the con- 
ventional framework of multiplicative renormalization for removing the ultraviolet 
divergence. Nevertheless, the SD equation for the ghost propagator is form-invariant 
under the multiplicative renormalization: 

We will return to the renormalization problem of the first horizon term later after 
discussing the second horizon term. 

3.1 The second horizon term 

It is well known that the composite operator/}^ [.c/]^ is multiplicative renormalizable 

D^,[^]^ = Z~^'\D^[s^]1^)r. (3.4) 
Hence the second horizon condition can be tamed by the conventional method: 

{h{x)) = - \im{{D^\s^]^)\D^\s^]^)'')k. (3.5) 

The unrenormalized horizon condition = (A^^ — 1)D — {h{0)) is 

= (iV2 _ ^^jj _ ^^(Q)^ _ _ - 1)[1 + n(0)]. (3.6) 

We can obtain the multiplicatively renormalized horizon condition as 

=Zc[{N' - l)D - (/i(O))] = (iV2 -1){D- l)Zc[l + n(0)] 
= {N'-1){D-1)[1 + ur{Q)1 (3.7) 

by adopting the renormalization relation: 

l+n(0) = Zc'[l+nR(0)]. (3.8) 

Therefore, we have the renormalized SD equation: 

GR\Q) = ^[l + URm+WR{Q), (3.9) 
8 



which is rewritten in favor of G{0) as 

Gr\0) = ^G-^'iO) - ^wniO)+WRiO), (3.10) 

where we have used the renormahzation: 

G-\0) = Zc'G],X0), w{0) = Zc'wr{0), (3.11) 
This renormahzation prescription is compatible with the renormahzed relation: 

G],'ie) = 1 + UR{e) + WR{e). (3.12) 

Thus the renormahzed horizon condition is satisfied only when 

ur{0) = -1, G],\0)=WRiO). (3.13) 

This is also obtained as a self-consistent solution of the SD equation. 

The scaling solution Gr{0) = oo follows immediately from wr{0) = which was 
assumed implicitly by previous studies. On the other hand, the decoupling solution 
Gr{0) < oo as a result of wr{0) ^ means the existence of the massless pole in the 
correlation function A^i,(/c). In fact, Wij(O) = is true at the tree level and holds in 
perturbation theory, since the massless pole coming from the elementary Faddeev- 
Popov ghost was already removed in the definition of the IPI function Xfj,u{k). For 
wr{0) 7^ to be realized, therefore, a massless boundstate must be formed anew in 
the channel by the non-perturbative effect and the massless pole must be generated. 
However, the existence of such a pole would not mean the appearance of the physical 
massless particle in the spectrum, since it is not a gauge-invariant object. Therefore, 
'Wr{0) 7^ does not contradict with the experiment. This consideration leads to a 
possibility of the decoupling solution satisfying the Kugo-Ojima criterion for color 
confinement. 

3.2 The first horizon term 

We return to the first horizon term 

{h{x)) = - lim((^?^^ X ^)^{g^^ X ^)^)fc, (3.14) 
fc— »o 

which leads to the unrenormalized horizon condition 
0={N^ -l)D-{h{0)) 
= {N^ - 1) {d[1 + n(0)] + w{0) - G(0)[n(0) + w{0)]^} . (3.15) 
However, this horizon condition is not multiplicatively renormalizable! 

=Zc[{N' - 1)D - {h{0))] 
= (iV2 - l)DZc + lim Zciig^^ x ^)^((?i4 x ^)^)fc 

= {N^ - 1) {dZc[1 + uiO)] + ZcwiO) - ZcGiOMO) + wiO)]^} 

^{N^ - 1) [d[1 + ur{Q)] + wr{Q) - G'^(0)[n^(0) + ^^.^(O)]^} . (3.16) 

Note that the renormahzation prescription 

I + u{Q) = Zc^[l + ur{Q)1 w{Q) = Zc^wr{Q), G'j,\Q)=Zc^G-\Q), (3.17) 

can not maintain the first horizon condition. 
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The conventional method: This is reasonable, since we did not use the localized 
renormalizable GZ theory. Without introducing the Zwanziger ghost field ^, ^, a;, a), 
the usual framework of the multiplicative renormalization does not work. Therefore, 
we must move to the localized GZ theory which is manifestly multiplicative renor- 
malizable. The Slavnov- Taylor identity means the horizon condition is translated to 
(See Appendix A) 

(n"'/V/^''"'<(x)er(x)) = {h{x)) = {N' - 1)D. (3.18) 

This is rewritten into the covariant derivative form, since the average does not depend 
on X due to the translational invariance: 

{tr'^'gD.l^f''^^^) = (MO)) = {N' - 1)D. (3.19) 

Then the horizon condition is multiplicatively renormalized: 

(z7;^'/'^?K(/^.H^^er)i?) = Zc{h{0)) = Zc{N' - l)D. (3.20) 

Thus an overall renormalization constant Zc is enough to renormalize this horizon 
condition. 

The effect of this horizon condition to the propagator of the Zwanziger ghost field 
uj^^{x) (see Appendix B) was calculated based on the localized Gribov-Zwanziger 
theory in section 10 of [22] • The Zwanziger ghost u has the same type of in- 
teraction to the gluon and hence it should have the same IR behavior as that of 
the Faddeev-Popov ghost. In this case, the actual calculation is very similar to 
that for the second horizon term. It was concluded that the propagator has a 
1/(A;^)^ singularity at k = 0. As far as the author understands, however, this con- 
clusion was derived with an implicit assumption (limfc^o[^^5'(^^)] = 0) similar to 
limfc^o[^^^(^^)] = linifc_>o w;(/c^) = w{0) = 0, i.e., absence of a massless pole in a cor- 
relation function (10.12) of [22] where g is assumed not to behave as gik"^) ~ c//c^ for 
/c^ ~ (c 7^ 0). If we remove the implicit assumption, therefore, we reach the same 
conclusion as that for the second horizon condition given in the previous subsection. 
Therefore, we do not pursue this direction. 

An unconventional method: Instead, we work in the original non-local Gribov- 
Zwanziger formulation. We recall the SD equation with the horizon condition: 

, , G^1(0)-2 + Ga(0) ^a(O) , , 

Gjy [k ) = — — ma(0) + UA{k ) + WA{k ). (3.21) 

If the horizon condition is incorporated into the SD equation, a partial cancellation 
at k = occurs between the horizon condition and the ghost self-energy. This cancel- 
lation always occurs for the multiplicative renormalizable part coming from A^^(O). 
This is not the case for the contribution from the remaining term A^^(O). 

In order to avoid the ultraviolet divergence, the UV cutoff A has been introduced 
thereby to make the self-energy part u{k'^) + w^k"^) finite. Consequently, G must 
depend on A, i.e., G{k'^,A) = G^ik'^)- In this sense, the above SD equation is an 
unrenormalized version. The multiplicative renormalization fails for the first horizon 
term which is non-linear in G~^, see (13.211) . 

The first horizon condition is not multiplicatively renormalizable. Therefore, the 
ultraviolet divergence can not be removed within this scheme. However, a novel 
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situation occurs by introducing the first horizon condition into the SD equation. By 
the resuhing combination ut^{k'^) — ma(0), the uhraviolet divergence cancels exactly 
and the ultraviolet cutoff A can be sent to infinity to obtain a finite function of /c^. 
Then the above SD equation is regarded as a self-consistent equation to give a finite 
ghost function G{k'^) = liniA^oo G'A(/i;^) < oo. The term w\{k'^) is finite from the 
beginning by some reason, although the A dependence is apparently assumed. See 
Appendix B. In other words, the SD equation is self-organized (in a non-perturbative 
way) to give a finite result. Therefore there is no need for the specific ultraviolet 
renormalization in the case of the first horizon term! Thus, the results obtained in 
the unrenormalized case hold also after the ultraviolet cutoff A is send to infinity, as 
far as the ghost propagator or dressing function is concerned. 

For the second horizon condition, a similar situation does occur. 

Gl\k'') = ^[l + uM] - «a(0) + u^{k^) + w^{k^), (3.22) 

which is also rewritten as 

G-^\e) = ^^(0) + u^^e) + w^ik'). (3.23) 

As already shown, due to the linearity, the UV renormalization of this SD equation 
can be performed within the multiplicative renormalization framework: 1 + ^^(0) = 
Zc^[l + ur{0)], wa{0) = Zc^wr{0) and Gl\0) = Zc^G^\0). Moreover, we can 
use the same argument as the above: The UV divergence cancels for m(/c^) due to 
the subtraction — ma(0) -|-ma(A;^), while wp^ik"^) is finite. Therefore, the UV cutoff in 
Gt^ik!^) is sent to infinity without divergences. 

Even if the ghost dressing function is ultraviolet finite, we have still finite renor- 
malization coming from the choice of the renormalization point. Therefore, we study 
the renormalization point dependence. From 

GR{k\^?) =Za\fi\A')G{k',A'), (3.24) 

we have 



Gnik^fi^) G{k^A' 



Gnifi^fi') G{^i\A^y 
In particular, at fc^ = and A = oo 



(3.25) 



G^(0,/i^) = G^,(/i^/i^)^%^ = gi?(/i^/i^) ^, , (3.26) 

G{fi'',oo) G[fi'',oo) 

When the decoupling solution with w{0) =0 is realized for the first horizon term, 

Gn{0,^') = Gn{^\^')^^ = G nif^' , f^')^^^ ■ (3.27) 

For instance, if one chooses the renormalization condition G Fi{fi^ , fi^) = 1 at = 
l.SGeV for D = 4, then G(/i^) = 1.2 at /x = l.SGeV for a given boundary condition 
G(0) = 3. In this way, we can reproduce the Orsay data [H]: Gij(0,/i^) = 3/1.2 = 2.5. 
The Sao Paulo data [12] with our interpretation G/j(0,yU.^) = 4.2 and G/j(/i^,yU^) = 2 
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at /i = IGeV is also consistent with this analysis G{^^) = 1.5 at yU = 1. Here the 
value G(yU^) is taken from the numerical solution of the Schwinger-Dyson equation 
or functional renormalization group equation with the boundary condition G{0) = 3 
given in [15] where it is claimed that it is only a matter of infrared boundary conditions 
G{0) whether scaling or decouphng occurs. 

4 Coulomb gauge 

We consider the Coulomb gauge 

djs^,{x) = (j = l,...,D-l), (4.1) 
where x = (x, x^) = (a^i, ■ ■ ■ x = (xi, ■ ■ ■ ,X£)-i)- Then the total action reads 

Otot c I C 

'-'yM • — '-'YM + '-'GF+FP, 

^YM := / d^'x^^^, ■ ^GF+FP := / d^'x [SS ■ dyS^j + i^i ■ djDj^} , (4.2) 
We can adopt the horizon term which is instantaneous in xd = yo '■= t: |27j 

/^[^(x)] = I rf^-V^?r^^<(x)(/^"^)^^(a:,Z/)^7r^^<(z/)U,=,,, (4.3) 



or 



h[^ix)] = I rf^-V^,K]^^(x)(i^-^)^^(x,y)D,K]^^(y)U,=,„ (4.4) 

where K is the Faddeev-Popov operator in the Coulomb gauge K[£/] := —djDj = 
—dj{dj + g^/jx). The horizon action Sh is non-local in space, but local in time. The 
horizon condition is given by 

{h{x))az = {N'-l){D-l). (4.5) 

Here {h{x))Qz is independent of t due to translational invariance. The horizon action 
Sh is non-local in space, but local in time. For instance, 

5^ =7 [ d^xh[£/{x)] =-f fdt f d^-^xh[£/{x)] 



=7 



JdtJ d^"ix J d^~'ygf^^^^fix)iK~'f^ix,y)gf^^^^fiy)\^,=y,. (4.6) 



The ghost propagator (i.e., the complete propagator of the Faddeev-Popov ghost) 
G^^{x — y) = ('^^(x)^'^(y)) is instantaneous, G'^^{x) = G^^{yi)6{t). In momentum 
space, it is independent of ko, i.e., 

^^^(k,^^) = G^^(k). (4.7) 

The free ghost propagator is independent of 

(k) = (4.8) 



12 



We can introduce the ghost dressing function (or form factor) (i(k^) which is dimen- 
sionless by 

G^^(k, ko) = G^^(k) = -5^^^. (4.9) 
In the similar way, we can show that the identity holds: 

rf(k2) = [1 + n(k2) + W7(k2)]-^ (4.10) 
It is also shown that the first horizon term (14.31) yields 
(/i(0)) = - lini((^i^j X ^)^{g£^j x ^)^)k 

= -{N^- 1) {{D - l)u{Q) + W7(0) - rf(0)[M(0) + w{Q)f] , (4.11) 

while the second horizon function (14.41) leads to 

(MO)) = - lim ((D,^)^(D,<^l^)k = -{N'- 1) {{D - 2)n(0) - 1} , (4.12) 



where M(k^) and w(k^) are defined by 
Af (k) := {{g^, X ^)\g^, x 



k2 



(4.13) 



Thus, we can obtain the similar results in the Coulomb gauge to those in the Landau 
gauge by replacing D in the Landau gauge with _D — 1 in the Coulomb gauge. We find 
one-parameter family of solutions including both the scaling and decoupling solutions. 
See e.g. [2H] for the ghost dressing function obtained by the numerical simulations on 
lattice. 

5 Conclusion and discussion 

In this paper we have discussed how the existence of the Gribov horizon affects the 
deep infrared behavior of the ghost propagator in the Landau and Coulomb gauge 
G = SU{N) Yang-Mills theory, using the Gribov-Zwanziger framework with the hori- 
zon condition {h{x)) = {dimG)D. Moreover, we have shown how to incorporate the 
horizon condition into the Schwinger-Dyson equation for the ghost propagator to 
discriminate between scaling and decoupling. We have examined two horizon con- 
ditions derived from two types of horizon terms, both of which were proposed by 
Zwanziger [20l[22]. We have shown that one parameter family of solutions parame- 
terized by w{0) exists in both cases, although some results crucially depend on the 
choice of the horizon term adopted. The value w{0) has been assumed implicitly to 
be zero w{0) = in previous studies. 
For the first horizon term 1201. 



h{x) = I d^ygf^^'^s^f{x){K-'f^ix,y)gf^^^^^^iy), (5.1) 

the GZ theory is not multiplicatively renormalizable. However, the SD equation for 
the ghost propagator and the ghost dressing function can be UV finite, once the 
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horizon condition is incorporated into the SD equation as proposed in this paper. 
The decouphng solution, i.e., finite ghost dressing function G{k'^) even in the hmit 
/c — *• 0: G'(O) < oo is allowed to exist, unless w{Q) = — oo. The Kugo-Ojima criterion 
m(0) = —1 is not necessarily satisfied except for a special choice of w{0) = 1/2 for 
any D, leading to G(0) = 2. For w{0) = 0, the boundary values are G{0) = 1 + VD 
and n(0) = {—D + \^) /{D — 1) up to renormalization point dependence. A possible 
renormalization scheme and the renormalization point dependence of the decoupling 
solution has been discussed. This should be compared with the paper [29] . 
For the second horizon term [22] 

h{x) = J d''yD,[s^f^{x){K-Y''{x,y)D,[s^f''{y), (5.2) 

the GZ theory is multiplicatively renormalizable. The Kugo-Ojima criterion u{0) = 
— 1 is satisfied both in the unrenormalized and the renormalized cases [301136] . For 
wr{0) = 0, the scaling solution, i.e., infinite ghost dressing function in the limit 
A; — > 0: Gji{0) = wji{0)'^ = oo, even after the renormalization. For Wfi{0) ^ 0, the 
decoupling solution is obtained against the claim in the previous literatures. 

Thus the investigation of wr(0) is crucial to see which solution is realized. An in- 
teresting step towards this direction was done in [29] . However, according to our anal- 
ysis, their result, namely, the decoupling solution with Wi?(0) = is not compatible 
with the multiplicative renormalization scheme they used: 1 -|-m(0) = Zq^[1 + ur{0)], 
w{0) = Zq^wr{0) and G~^{0) = Z^^G]i^{0). This issue should be reexamined and 
confirmed by further investigations. 

If we require that two horizon conditions give the same result, then the relation 
m(0) = —1/2 — w{0) must be satisfied for any D. This implies that for any D 

G(0) = 2, n(0) = -l, w{0) = ^. (5.3) 

Then, in the unrenormalized case, the Kugo-Ojima criterion is compatible with the 
decoupling solution for both choice of horizon term, against the conventional wisdom 
where the scaling solution is believed to be consistent with the Kugo-Ojima criterion. 
However, G(0), u{0) and w{0) are not renormalization group invariants. There is 
no guarantee that this relation is preserved. In fact, the Kugo-Ojima criterion is 
preserved only for the second horizon term. In order to judge which horizon function 
is realized, one need to know the result of numerical simulations on finer lattice, i.e., 
with smaller lattice spacing corresponding to larger ultraviolet cutoff, in addition to 
larger size lattices. But this might be impossible, because the horizon term alone 
does not uniquely fix the gauge, since there are still Gribov copies in the first Gribov 
region. 

Thus we can conclude that wr{0) plays the role of an additional non-perturbative 
gauge parameter which uniquely specifies the solution from a one-parameter family 
of solutions including the scaling and decoupling. In view of this, the value of wr{0) 
itself has no physical meaning. In fact, it has been shown [37] that all solutions 
(decoupling as well as scaling) lead to quark confinement by proving the vanishing 
of the Polyakov loop as a gauge-invariant order parameter of quark confinement. In 
this sense, discriminating between decoupling and scaling may not be so important 
from the physical point of view and main contribution to phenomenological studies 
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comes from the high- momentum region and the intermediate region around IGeV 
which is stable irrespective of adopting the scahng or decouphng solution. It is shown 
that both scaling and decoupling solutions do not contradict the general principles of 
quantum gauge field theories [38] . 

How the existence of the horizon is relevant for color confinement. In the Gribov- 
Zwanziger theory (restricted to the 1st Gribov region), the BRST symmetry is bro- 
ken by the existence of the horizon. = ^S-f ^ 0. Nevertheless, there exists a 
"BRST" like symmetry (without nilpotency [39] or with nilpotency [ID] which leaves 
the Gribov-Zwanziger action invariant. Then we could apply the Kugo-Ojima idea 
to the Gribov-Zwanziger theory, which opens the path to searching for the modi- 
fied color confinement criterion a la Kugo and Ojima. In view of this, defining a 
non-perturbative BRST transformation will be another interesting possibility to be 
investigated [H]. 

[Note added] 

In preparing this paper, it is pointed out in [13] that by writing u{0) and w{0) as 
function of G(0), ma(0) +oo and w\{0) — > — oo such that u\{0) + w\{0) — 1, 
provided that Ga(0) — *• oo as A — *• cxo, see Appendix C. This analysis uses the first 
horizon condition and the relation G{0)~^ = l + u(0) +w{0). However, the statement 
G'a(O) ^ oo as a ^ cxd is a result of perturbation theory. This analysis does not 
use the full information coming from the relation G(/c^)~^ = 1 + n(A;^) + w{k'^) (the 
Schwinger-Dyson equation for the ghost dressing function) for the whole momentum 
region. As we have shown in the text, Ga(0) remains finite even after the limit A — >■ oo 
in a non-perturbative way, once the full information is used. 
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A The function w and its ultraviolet behavior 



Two functions u and w are defined by 



which imphes 



1 



l)(iV2 
-1 



k^k" 



(D-l)(iV2_i) 



5^"^ -D 



A;2 

k^'k'' 



irAB 



k^ 



(A.l) 

(A.2a) 
(A.2b) 



The argument given in is based on the power counting. dim.[=g/] = [D — 2)/2 = 
dim.[^] = dim.[^] means dim.[^^i^^] = 2D - A. Amv.[g] = (4 - D)/2. Therefore, 
dim.[(^^i^^)j = 2D-A-D = D-A. Thus, for D = 4, \^l^{k) has at most 
logarithmic divergent. Only u has divergence, while w is ultraviolet finite. The IPI 
part has no massless pole. The ultraviolet divergence appears in u and not in w. 
So the ultraviolet divergence is proportional to g^p. In fact, the Brown- Pennington 



projector 



D 



fc2 



eliminates the term proportional to g^^. 



B A localized Gribov-Zwanziger theory 

The Gribov-Zwanziger theory can be rewritten into the local form [22] by introducing 
additional fields called the Zwanziger ghosts ^.^.uj.oj: 

^-,jd^.m ^ j[di][dl][duj][duj]e^v[-S,W , (B.l) 

where 

~S. =: / d^Ai'/K^^'f,' - -r^^^'^-r 

+ z7^/^^7/^^^<ef + ^7^/^^7/^^^<ef ]. (B.2) 
The localized action ^gz for the Gribov-Zwanziger theory is obtained 

=5ymM + 5gf+fpK, ^, m + ~S,W. I ^, cc], (B.3) 

where 

=5^GF+FP := / rf^'x [SS ■ d^s^^ + ■ d^D;^] . (B.4) 

The localized GZ theory is known to be multiplicatively renormalizable to all 
orders where the multiplicative renormalization factors are introduced as 

C/P ry^l'i'C/pR C/^ 171/2(^11 

To — Zjq To ^ To — Zjq To , 

g =ZggR, Zg = Z^Z'^I'^Z^^, (B.5) 
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and 



In the localized GZ action, the covariant derivative form can be used 



This means another form of the non-local horizon term 



(B.6) 

AC-^ACy 

(B.7) 
(B.8) 



h{x) = J d^yD[^]f{x){K-'f^{x,y)D[^]f{y) 

C u{0) and w{0) as functions of G'(O) 

The analysis done in |32] is peformed to two horizon terms to compare their implica- 
tions. 

For the first horizon term ^^^'^-'^ 



ua{0) 
wa{0) 



1 



- D-'{-{D - l)n(O) - G{0)[u{0) + w{0)]}, 
(MO)) 



D - 1 
1 



Ga(0) -1-D 



D-l ^ 
Using the horizon condition, 

ua{0) 
wa{0) 



Ga{0) + 2- D + D 



(MO)) 



Ga(0)' 



(C.l) 



1 



D- 1 

1 

D-1 



{Ga(0) 



D} 



{-Ga(0) + 2} + 



If Ga(0) — > 00 as a — >• cx), then ma(0) — >■ +00 and wa(0) 
wa{0) -1. 

For the second horizon term ^^'^^''^ 



(C.2) 

-00 such that Ma(0) -|- 



ua{0) 
wa{0) 



1 



^5P^-D-H-P-I)n(0) + 1} 
(MO)) 



D-1 
1 



1 - D 



D 



D-D 



D{N^ - 1) 
(MO)) 



Ga(0)- 



(C.3) 



Using the horizon condition. 



ua{0) 
wa{0) 



— 1 (A — indep. 
1 

Ga(0)' 



(C.4) 
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Figure 2: This graph is extracted from [43]. The solutions for n(0, A) and w;(0, A) plotted 
as a function of G(0,A). The particular solution, G{0) = 3, u{0) = —2/3, proposed by 
Kondo (black circles), obtained by imposing w{0,A) = 0, corresponds to the intersection of 
u + w and u. The current lattice solutions for the bare ghost dressing functions at vanishing 
momentum lie inside the green dotted square. The Kugo-Ojima parameter is in the region 
-1 < m(0) < 0, when 2 < G(0) < 5. 



If G'a(O) — s> oo as a — i> oo, then w\{0) such that ma(0) + w\{0) —1. It 
should be remarked that u\{0) = —1 independently of A. This result seems to be 
inconsistent with the lattice result [3n[32] . 

Remark: In the first horizon term, if the contribution from A^^ is neglected |l2], then 
j^^=D-^{-Du{0)-wm, 



ma(0) 
wa{0) 



1 



1 



D - 1 
1 



G'a(O) 



D 



(MO)) 



D - 1 

Using the horizon condition, 

ua{0) 
wa{0) 



-D + D 



D{N^ - 1 

(MO)) 



Ga{0) 



D 



1 



D - 1 
D 



GaIO) 



D <-l. 



D-IGa^QY 



(C.5) 



(C.6) 



If G'a(O) — > oo as a ^ oo, then ua(0) —1 and tfA(O) such that UA(0)+tfA(0) 
— 1. However, ma(0) < —1 which contradicts with the lattice result. 
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